On the quantum criticality in the ground and the thermal states of XX model 
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We compare the critical behavior of the ground state and the thermal state of the XX model. 
We analyze the full energy spectrum and the eigenstates to reconstruct the ground state and the 
thermally excited state. With the solutions, we discuss about several physical properties of the 
states, which are related to quantum phase transition, in various limits, at zero temperature as well 
as at a thermal equilibrium. 



PACS numbers: 

Conventional definition of quantum phase transition 
(QPT) is given by energy level crossingof ground energy 
state with varying external parameter [l| . The definition 
is based on the view that the ground energy state is taken 
as a quantum phase itself and, thus, QPT is a transition 
occurring at the point where a zero temperature state 
is transited from one state to the other. Because of the 
singular structure of the energy transition, the point of 
QPT is known to be detected by non-analyticity of free 
energy. 

Within the definition, a difficulty may arise if one tries 
to consider a finite temperature QPT. It is because QPT 
is only defined at T — 0. Typical treatment of finite 
temperature phase transition so far is provided by the 
non-analyticity of thermodynamic properties or the char- 
acteristic wave length of the system 0, However, it 
is not straightforward how the singularity in the thermal 
state is consistently combined with the notion of energy 
level crossing at zero temperature. 

In this article, we study the full solutions of XX spin 
chain model, not only the energy eignevalues but also the 
eignestates in the original spin basis. With the solutions, 
we demonstrate how the sequential energy level crossings 
at zero temperature is led to a continuous QPT in infinite 
spin limit, N — > oo. In the XX spin model, the number 
of ground energy level-crossings grow with the size of 
the chain. When N — > oo, the crossing points become 
dense within a sharply defined critical region, giving rise 
to a continuous QPT. This is a typical example when a 
ground state is no longer in a pure state but in a mixed 
state. At the same time, each crossing point of the ground 
energy state becomes fully analytic. 

The situation is comparable to a state in finite temper- 
ature which is a mixture of all the state spectrums with 
the Boltzmann distribution. We investigate the struc- 
ture of state in a thermally excited XX spin chain and 
the purity of the state to compare it with the ground 
state of the infinite chain. In the case, it can be seen 
that the temperature is acting as a controlling parame- 
ter to determine the degree of state mixture in the sys- 



tem. Thermodynamic properties of the thermal state are 
remained analytic in the region. As temperature is in- 
creased, the state moves into a thermal mixture which 
destroys the quantum correlations, entanglement, in the 
system. Extensive discussions on the relation between 
QPT and entanglement in the XX model can also be 
found at HHil. 

XX MODEL AND FULL SOLUTIONS 

We start our discussion by considering N spin 1/2 par- 
ticles in a line, coupled by nearest neighbor XX interac- 
tion, with Hamiltonian 

N 
i=l 

where we have assumed that the exchange interaction J 
fixes the energy unit. The model with periodic boundary 
condition has been firstly solved by Katsura 6] and re- 
cently discussed at finite size in Ref. Q ■ Here, we assume 
open boundaries, with <tn+i — 0. After Jordan- Wigner 
(JW) and Fourier transformations, 

' 1=1 \ 1 / m —l 

the Hamiltonian takes a diagonal form, H = 
Ek=i A kd{d k + NB where A fc = -2B + 
2 cos [(irk)/(N + 1)]. The 2 N eigenenergies are 
e i = Z)k=i ^k{ipi\d\d k \i>i) + NB and the corresponding 
eigenstates are found 

N 

M> = II(4) ai,) l n >. ( 3 ) 

k=l 

where the binary integer number a® € {0, 1} identifies 
eigenenergy as e\ = J2k=i Afc"/? + NB. The state |0) is 
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FIG. 1: (a) The spectrum of 2 4 energy levels for the N = 4 
XX spin chains. The energies vary as a function of mag- 
netic field B. Blue lines are the energy change of the ground 
state. Energy level crossings occur at the critical points in 
the ground state as well as in excited states. Ground state 
entanglement is jumping as the state is switched from an en- 
tangled state to another state, (b) The ground state energy 
per spin against the magnetic field B in the thermodynamic 
limit. The energy is plotted together with large spin chain 
case, N = 50. The continuous energy level crossings occur in 
the limit at the region < B < 1. 



rizlfc+i and its energy 
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as e" 



d k d k -i ■ ■ ■ di\ip°) 
„ - —(N - 2k)B - 2£f =1 cos (^) (refer 1 for 
detail). The index k is an integer in 1 < k < N and Bk 
represents the points where energy level crossings occur. 
When B < Bn or B > B\, no other crossing exists and 
the ground state is simply given by the fermion vacuum 
state. The ground state energy for N = 4 is plotted in 
Fig. [1] (a) by blue lines. 

The ground state in the spin bases after k 
level crossing can be obtained recursively \<Pg) — 



])® N where Sf 



^2/(N+ l)sin [(tvM)/(N + 1)] and the explicit form of 
the state is given as 



N+l 



E 



Chh—h\h,h) ■ ■ ■ , h) (4) 



h<l 2 <—<h 



the vacuum: dfc|f2) = OVfc. It is also useful to know that 
(ipi | if>j) = Si j. The energy diagram of e; for N = 4 as a 
function of B is illustrated in FigOJa). 



GROUND STATE OF XX SPIN CHAIN 

In this section, we discuss about the ground state of the 
XX model. The ground energy state is the lowest energy 
state which is a function of the transverse magnetic field 
B. Letting Bk = cos[kir/(N + 1)] and defining the fc-th 
region, B^+i < B < B^, we can find the ground state 



where |Zi, I2, •• • , h) is the state whose Zi-th, l%- 
thy • •, ifc-th indicate the position of flipped spin. 
The amplitudes in the state are Ci^...^ — 



rP(2) 



• sm 



P(k)h 



Here, the sum extends over all the permutation of the 
numbers from 1 to k, (denoted by P(k)). We point out 
that, at each crossing point (caused by the variation of 
B), the ground state jumps discontinuously from one 
symmetric subspace to the another, orthogonal to the 
previous one. To see the state explicitly, we illustrate 
the ground states for N = 4 in spin bases which are 



\<P° g ) 



- I T,T,T,T), \<P 1 g ) = <h\ l,T,T,T) + 4l T,l,T,T) + 4l U,U) + ar| T,T,T,I) 

= 02(1 i,i,r,t)+2| r,i,i,t) + 2| i,r,r,i) + i r,r,i,i)+V5| i,t,i,t)+V5| ra,ra; 

= 03I T, 1,1,4) + 41 i. t, i, I) + 41 4,1,1,4) + %! 4,4,4, T), 1^) = 1 4,4,4,4) 

I 



(5) 



where af 



1 ± 



the state, it is clear that the superscript of if indicates 
the number of down spins. The ground state at the fc-th 
region is composed of ArCfe = N\/[k\(N — k)l] orthogonal 
vectors in spin bases which are the possible choices of k 
spins, out of total N spins. The ground states are in an 
orthogonal subspace, (fglfg ) = 5k,k>- 



It is found that the ground energy spectrum and the 
ground state become continuous at the thermodynamic 
limit, N — > oo. In the limit, the sum in e k n turns into a 



From definite integral and thus the ground state energy reads 



lim ^?1= 2 - 

N^oo N 7T 



B (arccosB - ~) - \/l - B 2 



(6) 



Within the critical region the ground state energy is ana- 
lytic everywhere except for B = ±1. From the finite size 
analysis we observe that such a critical line becomes a 
dense set of level crossing (see Fig. Q] (b)) and therefore 
is possible to be considered as a line of continuous QPT. 
At the crossing points Bk, the state exists as a mixed 
state with equal mixture of the states in a neighboring 
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fc-th region, \<pg) and \<Pg +1 ), i.e. 



p{B k ) = -{\^){<p k g \ + \ Vg 



h+1 M +1 \) 



(7) 



Therefore, one can find that the state in the thermody- 
namics limit is also varying from a mixed state into the 
other mixed state continuously as a function of B. That 
is because the critical points is now continuous func- 
tion, B u = limjv^oo cos(kir/(N + 1)) = cos(7uj) where 
< oj < 1. Moreover, one would find that the mixedness 
of the state in the limit is remained as a constant for any 
B in -1 < B < 1 as Trp^) 2 = 1/2. The mixedness (or 
purity) 1/2 of the ground state in the limit can also be 
proved from the limiting case of a thermal state which 
will be discussed in the following section. 



Sr'.rSm'.m- The state |</?™) can be obtained from the l- 
th state \ipi) using the transformation in ([2]), |<^™) — 

r N 



t) 8Ar where 
(0 



y/2/{N + 1) sin [(irjk)/(N + 1)] and G {0,1}. The 

values of the label I for the binary vectors a; = J^ fe 4'^fc 
are determined by r and m through an index match 

/ = r + J^^Lo 1 nCs under the constraints J2k=i 4^ = m - 
(It starts from I = 1 when (r, m) = (1, 0).) Similarly, the 
probabilities of the Boltzmann distribution are given as 
pi = q" 1 — exp(— (3e" l )/Z. The lowest energy state for 
a fixed m become a ground state at a given region of 
B, i.e {(p™) = \<p g n )- Furthermore, the corresponding en- 
ergy eigenvalues by the subspace indexes m and r can be 
found as 



XX SPIN CHAIN IN THERMAL EQUILIBRIUM 

In this section, we investigate the state in a thermal 
equilibrium and discuss about the properties related with 
quantum phase transition. At a finite temperature, the 
state exists in a mixture of all the energy state distributed 
by Boltzmann statistics. The state in a diagonal basis 

= nf =1 (4) Q ^ I"), 4° e {°i U is s iven as 



(8) 



1=1 



where pi = exp(— j3e{)/Z is the Boltzmann distribu- 
tion of the thermal state with the partition function 

Z = e -(3NB (1 + e ~/3A fc ^ J n terms of the spin ba _ 

sis, the thermal state further can be decomposed into the 
states which are in a symmetric subspace as, 
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m— r— 1 



m I m\ / m 
Qr \<Pr IVPr 



(9) 



The index m in the superscript indicates a symmet- 
ric subspace where vector components of the state 
contain the m number of spins flipped and the sub- 
script index r identifies the ArC m = N\/[m\(N — m)l] 
orthogonal vectors within the subspace. Thus, the 
state satisfies the orthonormal condition, (ip" l \ip"> ) = 



N 



(N-2m)B-2 Vcos[7rfc/(7V+l)]4 r) - (10) 



fe=i 



This shows that the slop of the energy against B is in- 
variant for the states in a given symmetric subspace, with 
a fixed m, and, thus, it also gives a proof that the en- 
ergy level crossings of excited states occur only between 
the states in a different symmetric subspace. The level 
crossings of the excited states for N = 4 spins also can 
be seen in Fig. [1] (a). 
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FIG. 2: Distributions of thermal state population for the 
N = 2 XX spin chain when (a) kT = 0.1 and (b) kT = 1. The 
states ji/0 and \tp + ) are the singlet and the triplet. Critical 
behavior disappears as the system goes into a thermal equi- 
librium. 



As a special case, we inspect the density matrix of the 
thermally excited state when N — 2 which is in a form, 



J 



P?- 



= Pi\ T, TXT, T I + pi\ip~){i>~\ +P3|V' + )(V' + I +Pi\ U)U,I I 



(ii) 



r 



where |V' ± ) = ^75 (| t, I) ± | 4, ?))• The state is expended 
in the three different symmetric subspaces, where there 
is no spin flip (m=0), one spin flip (m=l) and two spin 



flips (m=2). Each subspace contains 2C0, 2C1 and 2C2 
orthogonal elements of states, {| 1\T)}j |V' + )} an d 

{| 4il)}- We plot the probability distribution of the 
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states pi in Fig[2j One would find that the distribution 
and the partition function become analytic for any B 
and T as T is away from zero. If we investigate entan- 
glement in the system, one would identify that the state 
become separable when ApiPA > (pi ~ P3) 2 , (from the 
negativity the negative eigenvalue of partially trans- 
posed state), leading to a critical condition for tempera- 
ture kT > 1.13459. Interestingly, the separability condi- 
tion for the two qubit thermal state is independent of the 
strength of magnetic field B except when T = 0. Thus, 
the singular behavior of entanglement with respect to B 
is washed out as the temperature increased. This also 
has been depicted by the purity of the two qubit system 
in Fig. H 



(b) 




d B Trp z 



B 



FIG. 3: (a) Purity of two qubit system as a function of mag- 
netic field B and temperature T. (b) Derivatives of the purity 
as a function of B at different temperatures, T = 0.01 (blue), 
T = 0.3 (pink), T = 0.8 (yellow) and T = 2 (green). At a 
low temperature region, the purity of the system shows peaks 
near the critical points B = ±1/2. This is the point where the 
derivative of purity w.r.t. B becomes singular when T = 0. 
At the high temperature region, the purity becomes mono- 
tonic convex function of B and the derivative of the purity 
approaches to a linear function. 

With the density matrix of the thermal state (jSJ), we 
investigate a physical property which may be directly re- 
lated to the critical behavior of the system. The mixed- 
ness of the thermal state for any N is 



Trp| 



m 

k=l 



1 - (1 + cosh/3A fc )- 



1. This proves that the ground state of the infinite XX 
chain in the region \B\ < 1 is the equal mixture of the 
neighboring states as (J7J). Furthermore, it is straightfor- 
ward to see that entanglement disappears in the region 
\B\ > 1 and at high temperature. The separability of in- 
finite XX chain with entanglement witness has been fully 
treated in Q. 

In our analysis of the XX spin chain, it has been 
shown that the fundamental difference between the ther- 
mal state and the ground state of the infinite chain is 
only the degree of state mixedness. In both of the cases, 
the changes of states are impossible to be detected by sin- 
gularity due to the continuity of state changes. In that 
sense, it can be pointed out that the energy level crossing 
is a particular mode of the ground state change when the 
state is in a pure state. In general, the state change can 
be detected by the change of entanglement, up to local 
unitary, which is destroyed by mixedness of the states. In 
fact, entanglement identifies the true quantum phase (or 
state) of a system since it is the property which can exist 
only in quantum systems whether they are pure or mixed 
[Io| . Our approaches are general enough to be used for 
the discussions about the QPT and related properties of 
any other many particle systems. 



Purity 




(12) 



FIG. 4: Plot of the purity as a function of B and T when 
N = 10. The state is in a pure state when Trp 2 = 1 and is 
in a completely mixed state when Trp 2 = 1/2 10 . The purity 
of state is singular at the critical points of B where the state 
exists as a mixture of two orthogonal symmetric states. It 
becomes smaller as T is increased. 



This characterizes how far the state is from a pure state. 
The purity is plotted in Figf5]for a small number of qubits 
(N = 2) and for a large number of qubits (N = 10) in 
Fig|4] At the T — > 00 limit {fi — > 0), the purity becomes 
\/2 N . Thus, in the high temperature region and the 
thermodynamic limit, N — > 00, the purity approaches 
to zero where the state is in a complete mixture of in- 
finitely many orthogonal states. At the opposite extreme, 
T — > 0, the purity becomes either 1 or 1/2 depending on 
the B value. It is 1/2 when B is cos[kir/(N + 1)] and 
1 otherwise. When N — > 00, the purity becomes 1/2, 
i.e. linijv^oo linix^o Trpy = 1/2 in —1 < B < 1 because 
there exits at least a single value of B in the interval to 
make A& = 3k. In the \B\ > 1 region, the purity is 
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